NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 


MISSILE  TERMINAL  GUIDANCE  AND  CONTROL 
AGAINST  EVASIVE  TARGETS 

by 

John  CS.  Swee 
March  2000 


Thesis  Advisor: 
Second  Reader: 


Robert  G.  Hutchins 
Harold  A.  Titus 


Approved  for  public  release;  distribution  is  unlimited 


STIC  QUALITY  INSPECTED  4 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing 
instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of 
information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden,  to  Washington  headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188) 
Washington  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

March  2000  Master’s  Thesis 

4.  TITLE  AND  SUBTITLE 

Missile  terminal  guidance  and  control  against  evasive  targets 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Swee,  John  CS. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  PERFORMING 

ORGANIZATION  REPORT 

NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  position  of  the 

Department  of  Defense  or  the  U.S.  Government. 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (maximum  200  words) 

Terminal  guidance  of  the  missile  against  an  evasive  target  is  explored.  The  two  main  types  of  guidance  laws 
employed  in  the  majority  of  missiles,  namely  proportional  navigation  and  command  to  line-of-sight  are  modeled  using 
Matlab0  Simulink™.  The  two-dimensioned  missile-target  intercept  geometry  is  simulated  for  a  point  mass  missile  and 
target.  Missile  velocity  due  to  its  drag  also  factored.  The  engagement  results  for  different  scenarios  with  target  doing  a  9- 
g  evasive  maneuver  are  then  compared  to  analyze  the  performance  of  hybrid  proportional  navigation  guidance  with  bang- 
bang  control. 

14.  SUBJECT  TERMS 

Missile  Guidance  Laws,  Proportional  Navigation,  Command  to  Line-of-Sight, 

15.  NUMBER  OF 
PAGES 

94 

16.  PRICE  CODE 

17.  SECURITY 
CLASSIFICATION  OF 
REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION  OF 
THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFI- CATION 

OF  ABSTRACT 

Unclassified 

20.  LIMITATION 

OF  ABSTRACT 

UL 

NSN  7540-01  -280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  239-18  298-102 


11 


MISSILE  TERMINAL  GUIDANCE  AND  CONTROL 
AGAINST  EVASIVE  TARGETS 


John  CS.  Swee 
Major,  Singapore  Navy 

B.Engg.,  Nanyang  Technological  University,  1994 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 

ELECTRICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
March  2000 


Department  of  Electrical  and  Computer  Engineering 


iii 


IV 


ABSTRACT 


The  ability  of  a  missile  to  intercept  a  target  in  its  flight  is  greatly  determined  by  the 
guidance  law  employed  in  the  guidance  processing  of  the  missile.  Two  main  types  of  guidance 
laws  are  employed  in  the  majority  of  missiles,  namely  proportional  navigation  (PN)  and 
command  to  line-of-sight  (CLOS).  The  effectiveness  of  CLOS  however  is  limited  to  shorter 
ranges  of  up  to  about  6km,  due  to  its  sensitivity  to  angular  tracking  errors  between  the  ground 
station  and  the  target.  PN  is  probably  the  most  widely  used  homing  guidance  law,  which  seeks  to 
null  the  line-of-sight  (LOS)  angle  rate  by  making  the  missile  turn  rate  be  directly  proportional  to 
the  LOS  rate.  PN  does  not  suffer  from  the  range  limitation  encountered  by  CLOS  because  it  is 
self-homing  and  relies  on  an  onboard  seeker  that  provides  target’s  LOS  information  directly.  We 
modeled  the  two-dimensioned  missile-target  intercept  geometry  with  CLOS  and  PN  guidance 
laws  using  Matlab®  Simulink™.  The  engagement  results  for  a  non-maneuvering  target  were  first 
established  as  a  benchmark  and  subsequently  compared  for  the  case  of  a  target  with  a  9-g  evasive 
maneuver.  While  conventional  PN  was  shown  to  be  effective  against  a  non-maneuvering  target, 
it  has  to  be  modified  to  improve  its  performance  against  a  maneuvering  target.  Simulations  for  a 
proportional  navigation  strategy  incorporating  bang-bang  control  were  carried  out  and  analyzed. 
The  performance  of  this  strategy  is  also  presented. 
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L  INTRODUCTION 


A.  BACKGROUND 

The  sequence  of  launching  a  missile  against  an  airborne  target  involves  various  stages. 
Of  particular  concern  for  this  study  is  the  terminal  guidance  of  the  missile  to  achieve  target 
interception.  The  various  methods  of  missile  guidance  can  be  broadly  categorized  under 
autonomous  (i.e.  self  or  homing)  guidance,  or  ground  commanded  guidance. 

Autonomous  Guidance.  This  type  of  guidance  requires  the  missile  to  carry  its  own 
guidance  system.  The  major  benefit  of  autonomous  guidance  is  the  ability  of  the  missile  to 
track  its  target  after  it  is  fired  and  frees  the  ground  control  station  to  perform  other  missions. 
This  feature  is  also  known  as  “fire  and  forget”.  The  missile  would  need  to  carry  a  seeker, 
which  invariably  makes  it  more  costly.  Generally,  proportional  navigation  is  employed. 

Ground  Commanded  Guidance.  With  ground  commanded  guidance,  the  computing 
power  resides  in  a  Fire  Control  System  (FCS)  at  the  launch  platform.  The  FCS  tracks  both  the 
target  and  missile  until  interception  occurs.  It  computes  the  target  trajectory,  determines  the 
necessary  missile  acceleration  and  transmits  the  guidance  commands  via  an  encoded  up-link 
to  the  missile.  This  makes  the  missile  less  costly,  but  it  is  range  limited,  as  the  tracking  error 
increases  with  target  range  from  the  FCS.  Generally,  ground  commanded  systems  employs 
command  to  line-of-sight  techniques. 

B.  OBJECTIVE 

The  ability  of  a  missile  to  intercept  a  target  in  its  flight  is  greatly  determined  by  the 
guidance  law  employed  in  the  guidance  processing  of  the  missile.  The  objective  of  this  study 
is  to  modify  the  guidance  law  to  improve  its  performance  against  a  maneuvering  target. 
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C.  RELATED  WORK 


It  is  known  that  proportional  navigation  (PN)  was  optimal  against  a  non-maneuvering 
target  [Ref  1  &  3],  but  was  not  effective  when  the  target  maneuvers.  Work  in  improving  the 
missile  guidance  laws  mostly  began  after  WWII.  The  closed  form  solution  for  PN  was  only 
recently  published  in  1990  [Ref  7].  Efforts  to  modify  the  PN  guidance  had  been  carried  out  in 
recent  decades,  and  new  algorithms  including  augmented  PN  have  been  proposed  to  improve 
its  performance  against  maneuvering  target  [Ref  10].  There  are  also  efforts  to  evaluate  the 
performance  of  missile  guidance  laws  against  a  maneuvering  target  [Ref  8,9,10  &  1 1],  but 
few  work  have  been  earned  out  with  Bang-bang  control  together  with  PN. 

D.  THESIS  ORGANIZATION 

For  this  study,  we  will  explore  the  performance  of  proportional  navigation  for  a  two- 
dimensioned  intercept  geometry  of  point  mass  target  and  missile. 

Chapter  II  provides  the  reader  with  a  basic  understanding  of  the  various  guidance  laws 
typically  employed  in  today’s  arsenal  of  missiles. 

Chapter  III  dwells  on  the  problem  formulation,  employing  proportional  navigation  for 
different  target-missile  intercept  scenarios  and  geometries.  We  first  assume  a  non- 
maneuvering  target  scenario  and  then  repeat  the  scenario  for  a  target  that  maneuvers  with  a  9- 
g  turn  about  2  seconds  prior  to  the  time  of  missile  interception. 

Chapter  IV  provides  the  details  on  the  simulation  model  used  to  evaluate  the 
performance  of  the  modified  proportional  navigation  guidance  proposed  in  this  thesis. 

Chapter  V  presents  the  results  and  analysis  of  the  simulations. 

Finally,  Chapter  VI  presents  the  conclusions  and  recommendations  of  this  study. 
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II.  MISSILE  GUIDANCE  LAWS 

A.  GENERAL 

This  chapter  explains  the  various  guidance  laws  that  are  typically  implemented  in  the 
missile  guidance  processing.  The  missile  guidance  system  provides  the  auto-pilot  (i.e.,  missile 
control  system)  with  the  necessary  lateral  acceleration  commands.  The  missile-target  intercept 
geometry  has  several  important  parameters  as  shown  Figure  1  below: 


Figure  1.  Missile-Target  Intercept  Geometry 


Typical  parameters  that  can  be  described  in  the  missile-target  intercept  geometry  are: 
Rm,Rt:  Tracker  to  missile  range.  Tracker  to  target  range  respectively 
R  :  Range  between  target  and  missile 

am,  at:  Tracker  to  missile  angle.  Tracker  to  target  angle  respectively 
a:  LOS  angle  between  missile  and  target 

Ym>  Yt:  Missile  velocity  vector  angle,  Target  velocity  vector  angle  respectively 
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The  following  paragraphs  will  describe  the  major  types  of  Command  Guidance  (i.e., 


Beam  Rider  and  Command  to  line-of-Sight)  and  Homing  Guidance  (Pure  Pursuit  and 
Proportional  Navigation). 


B.  BEAM  RIDER  GUIDANCE 

Beam  riding  guidance  is  one  of  the  simplest  form  of  command  guidance.  The  object  of 
beam  riding  (BR)  is  to  fly  the  missile  along  a  tracker  beam  that  is  continuously  pointed  at  the 
target.  A  typical  BR  geometry  is  shown  in  Figure  2  below. 


Figure  2.  Missile-Target  Intercept  Geometry 


From  Figure  2,  we  obtain  the  various  parameters  as  follows: 


(2.1a) 

(2.1b) 

(2.1c) 
(2.  Id) 
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The  crossing  range  error  (CRE)  i.e.,  the  distance  if  the  missile  from  the  beam  is: 

CRE  =  Rmsin(ot  -  otn)  .  (2.2) 

If  the  missile  is  always  on  the  beam  (i.e.,  CRE=0),  then  the  missile  will  surely  hit  the 
target.  Hence,  the  BR  guidance  law  to  drive  the  miss  distance  to  zero  is  to  make  the  missile 
acceleration  command  nc  proportional  to  the  CRE. 
nc  =  K.CRE 

=  K .Rms\n(at  -  ani)  .  (2.3) 

=  K.Rmipt  -  am)  for  small  (at  -  am) 

Hence,  we  observe  that  the  guidance  command  is  proportional  to  the  angular  error 
between  the  missile  position  and  the  tracker  beam.  However,  implementing  the  above  BR 
guidance  will  result  in  an  oscillatory  missile  acceleration  (Ref  [1]).  A  larger  value  of  K  (i.e. 
K=10)  will  give  a  smaller  miss-distance  but  the  missile  acceleration  oscillations  increases. 
Conversely,  a  smaller  value  of  K  (i.e.,  K=l)  will  have  less  oscillation  in  its  missile 
acceleration  but  a  larger  miss-distance. 

In  order  to  stabilize  the  BR  guidance  loop,  we  can  add  a  lead-lag  compensation 
network,  such  as  (Ref  [1]): 


G(s)  =  K  ’ 


1+5/2 

1+5/20 


=  10  K  ’ 


5  +  2 
5  +  20 


5  +  20 


C.  COMMAND  TO  LINE-OF-SIGHT  GUIDANCE 

Beam  riding  guidance  can  be  significantly  improved  by  taking  the  beam  motion  into 
account.  Adding  the  beam  acceleration  to  the  BR  guidance  equation  yields  command  to  line- 
of-sight. 
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Command  to  line-of-sight  (CLOS)  guidance  keeps  the  missile  in  the  LOS  between  the 
launch  point  and  the  target.  A  typical  flight  trajectory  of  a  CLOS  missile  is  shown  in  Fig.  3. 
The  distance  between  the  missile  and  the  desired  trajectory  (i.e.,  the  radar-to-target  LOS  line) 
is  defined  as  the  cross-range  error  (CRE).  The  ground  control  station  will  compute  and 
provide  necessary  missile  acceleration  commands  to  bring  the  CRE  to  zero. 


Figure  3.  CLOS  missile  flight  trajectory 


The  amount  of  error  in  the  CRE  at  the  point  of  intercept  is  dependent  on  the  range 
from  the  launch  point.  This  inherent  disadvantage  of  CLOS  restricts  its  use  to  short  ranges. 


The  angle  of  the  beam  with  respect  to  the  target  (see  Fig.  2)  is  given  by: 


at  =  tan  1  — 


(2.4a) 


Then  the  angular  velocity  and  acceleration  of  the  beam  are: 


Angular  Velocity  :  at 


d\  tan 


-1  y, 


*,y,  ~  y,x  r 
2  ,  2 
x,  +y, 

x,y,  -  y,*, 


(2.4b) 


Angular  Acceleration :  at 


where  R, 


at  cos  at-  a,  sin  at -2atRt 


^rjZty, 

R, 


and  a(>  are  acceleration  in  x  and  y  directions 


(2.4c) 


The  acceleration  perpendicular  to  the  beam  a,„,  can  be  expressed  in  terms  of  the 
inertial  coordinates  of  the  target  acceleration: 


-a,  sin  at  +  a.  cos  at 


=  /?,or  +  2i?;dt 


In  order  for  the  missile  to  stay  on  the  beam,  we  are  striving  to  ensure  that: 


at  and 


<r  -at 
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Hence,  the  commanded  missile  acceleration  perpendicular  to  the  beam  should  be: 
am„  =  Rmarn  +  2Rmarn 

=  Kot  +  2K<X  . (2.7) 

Adding  the  beam  acceleration  term  ,  amp  to  the  nominal  acceleration  generated  by  the 
beam  rider  equations  yields  the  CLOS  guidance.  Hence  for  CLOS,  the  missile  acceleration 
command  is  given  by: 

n.  =  K.CRE  +  am 

C  frlp 

=  K.Rmsm(at -am)  + Rmat +  2Rmat  ....  (2.8) 

=  K.Rm.(ot  -  am)  +  Rmat  +  2 Rmat  for  small  [at  -  am) 

D.  PURE  PURSUIT 

Pure  Pursuit  (PP)  guidance  seeks  to  keep  the  missile’s  velocity  vector  pointing  at  the 
target.  Typically,  it  has  relatively  simple  processing  avionics  and  works  well  only  against  very 
slow  targets.  Against  a  fast  moving  crossing  target,  the  lateral  acceleration  required  from  the 
missile  at  intercept  tends  to  be  very  high. 


Figure  4.  Pure  Pursuit  Geometry 


Figure  4  shows  PP  target-missile  intercept  geometry.  The  angular  error  between  the 
missile  vector  and  the  LOS  of  the  target  to  missile  (LOSmt)  is  shown  in  the  figure.  The 
guidance  algorithm  computes  the  missile  heading  to  bring  the  angular  error  to  zero.  PP 
guidance  is  not  effective  against  a  fast  target,  especially  when  the  target  maneuvers.  Hence,  it 
will  not  be  considered  in  this  study. 


E.  PROPORTIONAL  NAVIGATION 

While  CLOS  is  typically  employed  for  shorter  engagement  ranges;  for  longer  ranges. 
Proportional  Navigation  (PN)  guidance  is  preferred  and  some  form  of  seeker  (active  or 
passive)  is  built  into  the  missile  to  track  the  target. 


Figure  5.  Proportional  Navigation  Flight  Trajectory 
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PN  guidance  provides  acceleration  commands  to  the  missile,  which  are  proportional  to 
the  rate  of  change  of  the  LOStm,  i.e.: 

Commanded  acceleration,  nc  =  NYcdL 

The  PN  class  consists  mainly  of, two  kinds  of  guidance  laws  [Ref  8]: 

a)  For  True  Proportional  Navigation  (TPN),  nc  is  applied  normal  to  the  LOS. 

For  TPN,  the  effective  acceleration  component  normal  to  the  missile  vector  is 
given  as: 

nr 

=>  Commanded  acceleration  to  missile,  am  — . - .  .  (2.9) 

cos(0M  -eL) 

b)  For  Practical/Pure  Proportional  Navigation  (PPN),  nc  is  applied  normal  to  the 
missile  velocity  vector,  Vm. 

=>  Commanded  acceleration  to  missile,  am  =  nc  =  NYCQL  . (2-10) 

F.  BANG-BANG 

For  bang-bang  control,  instead  of  providing  an  acceleration  command  based  on  a 
proportionality  relationship,  the  missile  maneuvers  by  ‘banging’  to  its  designed  g-limit  to 
bring  the  missile  along  the  LOS. 

Thus  the  commanded  acceleration  based  on  the  ‘bang-bang’  algorithm  is: 

Commanded  acceleration,  nc  =  glimjt  *  9.81  *  sign(Vc0L )  . (2- 1 1) 
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III.  PROBLEM  FORMULATION 


A.  MISSILE  TARGET  SCENARIO  AND  GEOMETRY 

For  this  study,  we  examine  the  performance  of  proportional  navigation  for  a  two- 
dimensioned  intercept  geometry  of  point  mass  target  and  missile.  The  target  and  missile  are 
assumed  to  be  point  mass  models  in  a  plane,  moving  with  velocities  VT  and  Vm  respectively. 

The  initial  position  of  the  missile  is  assumed  to  be  the  reference  point  of  the  relative 
coordinate  system  with  its  initial  velocity  vector  pointing  at  the  initial  target  position. 


Figure  6.  Proportional  Navigation  Geometry 


From  Figure  6,  we  obtain  the  various  parameters  as  follows: 
Angular  components: 
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dM  =  tan" 


0T  -  tan _1  — 
x. 


(3.1b) 


(3.1c) 


Closing  Velocity: 


dR  _ 


-VT  cos (0T  -dL)  +  VM  cos(0M  -  eL ) 


(3. Id) 


Alternatively,  using  pythagoras  theorem: 


R2  =  Ax2  +A y2 


2(Ax)Ax  +  2(Ay)Ay 

(Ax)Ax  +  (Ay)Ay 
— — 


Since 


we  have: 


A x  =  R  cos(0L ) 
A y-  R  sin(0L  ) 


V  - =  -[Ax  cos(0,  )  +  Ay  sin(0L )] 


(3.1e) 


LOS  rate: 


d\  tan 


-A  yt-y„ 


x,  ~  xm 

t  m 


dt  dt 

Ay  cos(0L )  -  Ax  sin((9L ) 
~  R 


(3. If) 
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B.  TARGET  AND  MISSILE  MANUEVER 


Turn  rate  of  an  aircraft  is  a  function  of  it’s  speed  and  acceleration  (Ref  [1]).  The 
maneuver  by  target  and  missile  in  this  study  is  implemented  by  using  turn  rates.  Fig.  6  shows 
an  aircraft  moving  at  a  constant  speed  along  a  circle  of  radius  R. 


Aircraft 
(point  mass) 


Figure  7.  Aircraft  (point  mass)  during  a  turn 


From  this  we  obtained  the  continuous-time  state  equation  as: 


x(0  = 


x(t) 

'0 

1 

0 

0 

Vx(t) 

0 

0 

0 

-UJ 

m 

0 

0 

0 

1 

Vy 

0 

m 

0 

0 

=  A(m)x(t) 


x(r) 


(3.4) 


C.  GUIDANCE  LAW  IMPLEMENTATION 


In  this  study,  we  study  the  performance  of  proportional  navigation  with  N’=3  in 
combination  with  bang-bang  guidance,  subject  to  a  20-g  limit. 

Proportional  Navigation:  The  proportional  navigation  guidance  Law  is  given  by: 
Commanded  acceleration,  nc  =  N’VC6L 

a)  For  True  Proportional  Navigation  (TPN),  nc  is  applied  normal  to  the  LOS. 


Commanded  acceleration  to  missile,  am  =  - 


b) 


cos(eM  -eL) 

For  Practical/Pure  Proportional  Navigation  (PPN),  nc  is  applied  normal 


to  the  missile  velocity  vector,  Vm. 


Commanded  acceleration  to  missile,  a 


We  know  that  turn  rate,  m  =  - 


acceleration 


perpendicular 


||Velocity|| 


=>E7„ 


=  nc 


Thus  the  corresponding  commanded  turn-rate  is  given  by: 


a) 


For  TPN: 


- 


vM  cos(0M  -  eL) 


b)  For  PPN: 


NjA 

'  ^ 


(3.5) 


(3.6) 
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where: 


N’  is  the  proportional  constant,  which  typically  varies  from  3-5. 

Vc  is  the  closing  velocity  obtained  from  Eqn  3. 1(d)  or  3. 1(e) 

0  L  is  LOS  rate  obtained  from  Eqn  3.1(f) 

Vm  is  the  instantaneous  missile  speed 

Bang-Bang:  For  Bang-bang  control,  instead  of  a  proportionality  relationship,  a  full 
20-g  acceleration  in  either  direction  of  the  LOS  is  applied  to  null  the  LOS  rate  [Ref  6].  This 
is  implemented  by:  Commanded  acceleration,  nc  =  20  *  9.81  *  sign(Vc0L ) 

Guidance  Laws’  Implementation.  This  study  examines  the  following  schemes  of 
guidance  laws: 

a)  TPN  guidance  alone 

b)  Bang-bang  guidance  alone 

c)  Hybrid  TPN  and  Bang-bang  guidance 

i)  Start  with  Bang-bang  and  switches  to  TPN  at  R<2km 

ii)  Start  with  TPN  and  switches  to  Bang-bang  at  R<2km 


D.  TURN-RATE  TIME  CONSTANTS 


A  parameter  associated  with  missile  (as  well  as  target)  maneuverability  is  the  turn-rate 
time  constant,  TWm-  The  relationship  of  the  turn-rate  time  constant  is  given  as  : 


GJ 


output 


=^1+T-S) 

=mc  (l+T^s) 


(3.7) 
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Generally,  T,Um  increases  with  altitude  and  decreases  with  missile  velocity.  A  fast  turn- 
rate  time  constant  with  a  large  navigation  constant  may  result  in  instability  in  the  overall 
guidance  system  [Ref  1].  Typical  time  constant  varies  from  0.5- 1.0  seconds. 

E.  MISSILE  DRAG 

The  missile  speed,  Vm,  will  encounter  atmospheric  frictional  drag  during  its  flight. 

This  will  decrease  its  speed  during  the  coast  phase  of  its  flight.  The  atmospheric  drag  is 
computed  as  follows  [Ref  1]: 

Drag  =  p*CD*  Area*  VM2/2  .  (3.8) 

where: 

Vm  =  Missile  speed 
Area  =  Reference  area 

(Which  can  be  taken  to  be  7t*radius2  for  a  cruciform  missile  body. 

We  assume  a  medium-sized  300  kg  missile  with  a  radius=0.15m) 

Cd  =  Coefficient  of  Drag 

(Cd  =  0.2  was  chosen  for  supersonic  speed  during  its  coasting  flight) 
p  =  Atmospheric  density 

i)  For  height  <  9144m,  p=1.22557*exp(-h/9144) 

ii)  For  height  >  9144m,  p=1.75228763*exp(-h/6705.6) 

Missile  engagement  against  an  aircraft  typically  occurs  below  9144m.  For  this 
simulation,  an  arbitrary  height  of  2000m  was  used. 

The  drag  expression  given  in  Equation  3.8  applies  to  a  frontal  atmospheric  drag.  We 
can  expect  the  drag  to  increase  as  the  missile  does  a  turn,  since  the  effective  area  encountering 


atmospheric  friction  in  the  direction  of  the  missile  velocity  vector  increases.  The  coefficient 
of  drag,  Cd,  increases  as  a  function  of  coefficient  of  lift  Cl.  The  relationship  is  (Ref  [2]  &  [9]): 

Cd  =  Cdo  (1+ki  Cl2) 

where  : 

Cdo  is  the  nominal  coefficient  of  drag  with  angle  of  attack  is  zero. 

ki  =  a  constant  unique  to  the  aerodynamic  surfaces  of  the  missile 

Cl  =  coefficient  of  lift  for  corresponding  angle  of  attack 

We  attempt  to  establish  a  relationship  of  tum-rate  on  the  missile  drag.  We  assume  that 
when  the  missile  is  doing  its  maximum  tum-rate  at  20-g,  the  drag  coefficient,  Cd,  increases  by 
a  factor  of  5.  Following  the  above  equation,  a  similar  relationship  of  Cd  due  to  missile  turn 
rate  can  be  assumed  to  be  described  by: 

CD  =  CDo(l+k2o2)  .  (3.9) 

where : 

cr  =  missile  tum-rate 

k2  =  4/tn2 ^  ,  assuming  that  at  maximum  tum-rate,  Cd  =  5Cdo 

F.  MISSILE  VELOCITY  COMPENSATION 

Since  the  missile  is  expected  to  encounter  deceleration,  which  causes  the  missile  speed 
Vm  to  change,  it  is  reasonable  that  some  improvement  might  be  gained  by  compensating  for 
the  velocity  change  (Ref  [4]).  The  conventional  TPN  guidance  is  modified  as  follows: 

a„=N'V A-V.sin^-^)  .  (3.10) 
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G. 


SCENARIOS  OF  MISSILE-TARGET  INTERCEPTION 


The  scenario  adopted  in  this  study  assumes  an  inbound  target  flying  about  500m/s  (i.e. 
about  Mach  1.5).  The  missile  is  assumed  to  coast  at  an  initial  speed  that  is  twice  the  target’s 
speed,  i.e.  lOOOm/s.  The  initial  positions  of  the  missile  and  target  were  selected  so  that  the 
interception  range  is  about  6-6. 5km. 

The  following  scenarios  used  to  evaluate  the  performance  of  the  Guidance  scheme  are: 
Scenario  #1:  Crossing  target  with  lateral  distance  of  1000m  (See  Fig.  8). 

Part  A:  Target  on  straight  course 

Part  B:  Target  does  a  9-g  turn  away  from  missile  at  about  2  seconds 
prior  to  interception 


Cross-Range  Scenario 


x-coordinates  (meters) 

Figure  8.  Scenario  #1:  Crossing  Target  with  lateral  distance  of  1km 
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Scenario  #2:  Crossing  Target  at  45°.  (See  Fig.  9). 


Part  A:  Target  on  straight  course 

Part  B:  Target  does  a  9-g  turn  away  from  missile  at  about  2  seconds 
prior  to  interception 


Figure  9.  Scenario  #2:  45°  Crossing  Target 


Target  Initialization: 


Scenario#l: 

Lateral  Cross-Range 

Scenario#2: 

45°  Cross-Range 

Xt 

10000 

8000 

Vx, 

-500 

-500sin(45°) 

yt 

1000 

800 

Vy, 

0 

500sin(45°) 
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Missile  Initialization: 


Scenario#! : 

Lateral  Cross-Range 

Scenario#2: 

45°  Cross-Range 

0 

0 

Vxm 

200 

200 

ym 

0 

0 

Vym 

500 

500 

H.  NOISE  IN  LOS  RATE 

The  key  parameter  in  the  proportional  navigation  guidance  law  is  the  LOS  rate.  The 
performance  of  the  different  guidance  schemes  is  first  evaluated  without  noise  in  the  system. 
Subsequently,  a  sensitivity  analysis  was  carried  out  to  evaluate  the  performance  of  the 
proposed  guidance  strategy  by  introducing  an  additive  white  gaussian  noise  in  the  LOS  rate. 

The  performance  of  the  guidance  strategy  was  evaluated  for  a  range  of  standard 
deviations  in  the  noisy  LOS  rate  from  0.001  to  0.05  radians/sec  (0.06  to  2.86  deg/s). 


20 


rv.  SIMULATION 


A.  OVERVIEW  OF  MODEL 

The  guidance  system  simulation  model  is  generated  using  MATLAB®  and 
SIMULINK™.  It  comprises  of  the  following  blocks: 

1 .  Initialization 

2.  Target  Dynamics 

3.  Missile  Dynamics 

4.  Guidance 

The  above  simulation  blocks  are  described  in  the  following  paragraphs  in  this  Chapter. 
The  outputs  of  the  simulation  are  target  and  missile  state  vectors,  LOS  rate,  and  missile 
acceleration  commands. 

B.  INITIALIZATION 

Before  the  simulation  can  be  run,  the  initial  target  and  missile  state  vectors,  as  well  as 
the  sampling  time  and  duration  of  simulation  must  be  known.  This  initialization  is  given  in 
the  Matlab  file  ‘Init.m’.  The  following  variables/vectors  are  initialized: 

•  tinit :  target  state  vector  for  the  respective  scenario 

•  minit :  missile  state  vector  for  the  respective  scenario 

•  sampling  time 

•  Drag  component,  beta  =  p*Co*Area*  /2 

•  V]aunch'  closing  speed  measured  at  launch 
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c. 


TARGET  DYNAMICS  BLOCK 


The  tum-rate  tu,  is  given  as  a  step  function  with  an  amplitude  equal  to  9-g  (i.e.,  9*9.81 
m/s2).  This  is  used  to  give  the  target  a  9-g  evasive  maneuver  at  about  2  seconds  from  the 
expected  time  of  intercept.  It  has  a  tum-rate  time  constant  Tt=  1  second. 

From  Equation  3.4,  the  target  state  equation  is: 

'0  1  0  O' 

0  0  0  — <57 

=  0  o  0  1  x,arg£'° 

0  xu  0  0 

t  arg  et  J 

=  A(m)xtaiget(t) 

This  is  performed  by  the  Matlab®  function  ‘Targetdynamics.m’,  which  receives  w  and 
x(r)  as  inputs,  and  provides  x(t)  based  on  Equation  3.4  as  its  output. 

The  integrator  block  is  initialized  to  the  target  initial  state  vector  ‘tinit’. 


D.  MISSILE  DYNAMICS  BLOCK 


Figure  11.  Missile  Dynamics  Block 


The  missile  dynamics  block  is  similar  to  the  target  dynamics  block  with  a  tum-rate 


time  constant  Tm  =  0.25  to  1.0  seconds 

The  missile  state  vector  is  coded  in  the  file  ‘Missiledynamics.m’,  which  also  includes 
the  missile  drag  component. 

Missile  Drag.  The  missile  drag  component  is  implemented  as  follows: 

a) .  Atmospheric  drag  is  given  as: 

Drag  =  beta  *Vm2 

where 

beta=p*Co*Ref_area  /2  (initialized  in  Tnit.m’) 

Since  Force=mass*acceleration,  the  deceleration  due  to  drag,  adms  is: 

admS  =  Drag/masSmissiie 

=  beta  *Vm2/massmissiie 
The  corresponding  tum-rate  due  to  adrag  is: 

G?drag  =  &drag  /Vm 

=  beta  *  V^/masSmissiie 

b) .  Additional  drag  due  to  missile  tum-rate  tu,  is  assumed  to  have  a  relationship: 

k2  OJ  &Fdrag 

where 

k2  =  4/tn2max,  assuming  that  at  maximum  tum-rate,  the  total  drag 
component  increases  to  5  times  the  nominal  drag. 


23 


E.  GUIDANCE  BLOCK 


Additive  White  Noise 


INPUT:  Target  State 
(From  Target  Dynamics) 


w:  OUTPUT 

(To  Missile  Dynamics) 


(From  Missile  Dynamics) 


Figure  12.  Missile  Dynamics  Block 


The  missile  guidance  block  comprises  the  following: 

1)  Inputs:  The  input  to  the  guidance  block  are  the  target  and  missile  state  vectors 

2)  Proportional  navigation  variables:  The  intercept  geometry  is  obtained  from  the 
target  and  missile  vectors.  The  Matlab®  file  ‘PNvariables.m’  calculates  the 
respective  variables  used  in  proportional  navigation  as  described  in  Equations 
3.1(a) -3.1(e). 

3)  Guidance  Algorithm:  The  guidance  law  is  coded  in  this  block.  The  different 
guidance  strategies  are  coded  in  the  following  Matlab®  files: 

a.  ‘PNGuidance.m’:  Conventional  proportional  navigation  with  N’=3 

b.  ‘PNBangbang.m’:  20-g  bang-bang  guidance 

c.  ‘BangPN2km.m’:  Starts  with  bang-bang  and  switches  to  proportional 
navigation  when  the  range  of  target-to-missile  is  less  than  2km 
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d.  ‘PNbang2km.m’:  Starts  with  Proportional  navigation  and  switches  to 
bang-bang  when  the  range  of  target-to-missile  is  less  than  2km 
4)  Velocity  Compensation:  An  additional  turn  rate  to  compensate  for  the  missile’s 
velocity  changes  is  coded  in  the  Matlab®  file  ‘Velcomp.m’. 
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V.  PERFORMANCE  EVALUATION 


A.  PROPORTIONAL  NAVIGATION 

In  an  ideal  case,  where  the  missile  is  assumed  to  respond  instantaneously  to  the 
guidance  commands  (i.e.,  turn-rate  time  constant,  Tm  =T,  =  0),  conventional  proportional 
navigation  guidance  (PNG)  performs  optimally  for  a  non-maneuvering  target.  In  practice,  the 
missile,  as  well  as  the  target,  are  expected  to  have  turn-rate  time  constants.  In  addition,  the 
missile  which  is  not  in  sustained  flight  (i.e.  propulsion  is  burnt),  will  have  its  speed  reduced 
by  aerodynamic  drag. 

However  when  turn-rate  time  constants  are  considered  in  the  guidance  system,  PNG 
resulted  in  significant  miss  distances  when  the  time  constant  is  large  and  especially  when  the 
target  does  a  9-g  turn  away  from  the  missile  during  the  terminal  phase.  We  assumed  a  typical 
tum-rate  time  constant  Tt  =  1  second,  for  the  target  and  varied  the  missile  turn-rate  time 
constant,  Tm,  between  0.25  and  1  seconds. 

The  simulation  was  run  for  missile  drag  due  to  atmospheric  friction  alone  (Tables 
1(a))  and  also  for  the  case  of  additional  drag  due  to  the  missile  tum-rate  (Table  1(b)).  Detailed 
results  of  the  simulations  are  given  in  Appendix  C.  A  summary  of  the  miss  distances  is 
tabulated  below  for  the  two  scenarios  considering  different  missile  tum-rate  time  constants. 


Time  constant, 
Tm  (sec) 

Miss  Disl 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

0.25 

0.5096 

6.0796 

0.3271 

0.1018 

0.5 

18.6431 

0.2782 

3.0831 

0.6741 

34.1685 

6.8659 

4.3064 

Table  la.  Miss  Distances  for  PNG  (Tt  =  1.0s  and  Tm  =  0.25s,  0.5s,  1.0s)  without 
considering  drag  due  to  missile  tum-rate 
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Time  constant, 
Tm  (sec) 

Miss  Disl 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

0.25 

0.6468 

7.5479 

0.3601 

1.2425 

im 

0.1290 

19.8792 

0.1381 

4.5232 

| 

0.6207 

34.7688 

2.6993 

5.3368 

Table  lb.  Miss  Distances  for  PNG  (Tt  =  1.0s  and  Tm  =  0.25s,  0.5s,  1.0s)  with 
additional  drag  due  to  missile  turn-rate 


From  the  results,  we  observed  that  the  missile  time  constant  has  a  significant  impact 
on  the  miss  distance.  A  larger  time  constant  makes  the  missile  sluggish  and  less  able  to 
respond  to  the  given  guidance  commands  especially  when  the  target  maneuvers. 

Scenario#l  presents  a  more  difficult  target  to  intercept  than  scenario#2,  resulting  in 
larger  miss  distances.  This  is  because  the  geometry  of  scenario#l  requires  the  missile  to 
perform  a  larger  lateral  acceleration  during  the  terminal  flight  as  it  approaches  the  target. 

We  observed  that  Tm  =  1.0s  gave  unacceptably  large  miss  distances.  Thus,  the  missile 
time  constant  should  be  less  than  1  second  for  the  above  scenarios.  For  the  subsequent 
simulations,  we  assume  Tt  =  1.0s  and  Tm  =  0.5s  and  explore  modifications  to  the  guidance 
strategy  to  seek  improvements. 


B.  VELOCITY  COMPENSATED  PROPORTIONAL  NAVIGATION 

The  aerodynamic  drag  simulated  in  our  model  caused  a  reduction  in  missile  velocity 
of  about  150  m/s  over  7  seconds  of  missile  flight  (see  Fig.  11a),  when  we  only  consider 
frontal  atmospheric  drag. 

With  the  additional  drag  due  to  missile  turn-rate,  the  missile  velocity  profile  is  similar 
to  Fig  11a  when  the  target  does  not  maneuver.  However,  the  missile  velocity  reduces  sharply 
when  the  missile  turn-rate  increases  to  respond  to  target  maneuvers  (see  Fig  1  lb). 
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Missile  vector  velocity,  Vm  (meters/sec 


time  (seconds) 


Figure  11a.  Missile  velocity  for  scenario#!  without  considering  drag  due  to  turn-rate 


Fig.  3  Missile  Velocity 


Figure  lib.  Missile  velocity  for  scenario#!  taking  into  consideration  drag  due  to  turn-rate 


Velocity  compensated  proportional  navigation  (VCPN)  is  a  modification  of  PNG  to 


compensate  for  the  change  in  missile  velocity  by  an  additional  acceleration  command 


K  sin(0M  -0L). 

Detailed  results  of  the  all  the  simulation  runs  in  this  study  are  given  in  Appendix  C. 

A  summary  of  the  miss  distances  is  tabulated  below  for  the  two  scenarios  considering 
different  drag  model. 


Miss  Disl 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

0.4969 

18.6431 

0.2782 

3.0831 

VCPN 

0.5420 

18.1240 

0.0605 

Table  2a.  Miss  Distances  for  VCPN  guidance  (Tt  =  1.0s  and  Tm  =  0.5s)  without 
considering  drag  due  to  missile  tum-rate 


Miss  Disl 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

0.1290 

19.8792 

0.1381 

4.5232 

VCPN 

0.0480 

19.0654 

0.1240 

3.9718 

Table  2b.  Miss  Distances  for  VCPN  guidance  (Tt  =  1.0s  and  Tm  =  0.5s)  with 
additional  drag  due  to  missile  tum-rate 


From  the  results,  VCPN  guidance  does  not  offer  much  improvement.  This  is  not 
surprising  as  the  velocity  compensation  component,  Vm  sin(0M  small,  especially 

when  the  missile  lead  angle  is  often  less  than  a  few  degrees.  Furthermore,  since  any  change  in 
the  missile  velocity  will  be  reflected  in  the  closing  velocity,  Vc,  the  conventional  PNG  which 
computes  its  acceleration  commands  using  Vc,  in  essence  will  address  the  change  in  missile 
velocity. 
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C.  BANG-BANG  GUIDANCE 

For  PNG,  the  responsiveness  of  the  guidance  against  target  maneuvers  is  improved 
with  larger  proportional  navigation  constant.  This  led  to  considering  a  bang-bang  guidance 
strategy. 

Detailed  results  of  the  simulations  are  given  in  Appendix  C.  A  summary  of  the  miss 
distances  is  tabulated  below  for  the  two  scenarios  considering  different  drag  models. 


Miss  Disl 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

0.4969 

VCPN 

0.5420 

18.1240 

2.8216 

Bang-bang 

0.1116 

0.3318 

0.5761 

0.5214 

Table  3a.  Miss  Distances  for  Bang-bang  guidance  (Tt  =  1.0s  and  Tm  =  0.5s)  without 
considering  drag  due  to  missile  turn-rate 


Miss  Disl 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

0.1290 

19.8792 

0.1381 

4.5232 

VCPN 

0.0480 

19.0654 

0.1240 

3.9718 

Bang-bang 

0.4433 

0.5531 

0.0359 

0.0854 

Table  3b.  Miss  Distances  for  PNG  (Tt  =  1.0s  and  Tm  =  0.5s)  with  additional  drag  due 
to  missile  turn-rate 


The  simulation  results  showed  that  the  performance  of  Bang-bang  guidance  against  a 
maneuvering  target  was  much  better  than  both  conventional  PNG  and  VCPN. 

However,  bang-bang  is  sensitive  to  changes  in  sign  of  LOS  rate,  which  occurs  when 
the  guidance  law  tries  to  null  the  LOS  rate.  This  makes  bang-bang  inherently  more  susceptible 
to  noisy  LOS  rate. 
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In  scenario#l,  which  presents  more  frequent  sign  changes  in  LOS  rate  about  zero,  we 
observed  that  the  amount  of  control  effort  used  by  bang-bang  was  considerably  higher  at 
about  144-149  compared  to  75-76  for  the  case  of  PNG  and  VCPN  (see  Table  4). 


Missile  Control  Effort  (Integration  of  turn-rate) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

76.2543 

200.9011 

622.9764 

711.1141 

VCPN 

75.502 

199.255 

608.6091 

681.0754 

Bang-bang 

149.3850 

243.9412 

671.2259 

701.4282 

Table  4a.  Missile  Control  Effort  (Tt  =  1.0s  and  Tm  =  0.5s)  without  considering  drag 
due  to  missile  turn-rate 


Missile  Control  Effort  (Integration  of  turn-rate) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

76.2628 

196.771 

651.076 

736.5428 

VCPN 

75.4792 

196.0564 

609.9731 

698.861 

Bang-bang 

144.6337 

240.7891 

655.8412 

707.0795 

Table  4b.  Missile  Control  Effort  (Tt  =  1.0s  and  Tm  =  0.5s)  with  additional  drag  due  to 
missile  turn-rate 


Bang-bang  has  a  disadvantage  of  frequent  swinging  of  acceleration  commands  when 
LOS  rate  is  near  zero,  requiring  larger  missile  control  efforts. 


D.  COMBINATION  OF  PROPORTIONAL  NAVIGATION  WITH  BANG-BANG 

It  is  conceivable  that  we  can  employ  the  advantage  of  proportional  navigation 
guidance  with  Bang-bang.  Hence  a  combined  strategy  of  VCPN  and  bang-bang  guidance  was 
explored.  Two  strategies  were  examined: 

a.  Bang_PN:  Begins  with  bang-bang  and  switches  to  proportional  navigation 
when  the  missile  approaches  within  2km  from  the  target. 
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b.  PN_Bang:  Begin  with  proportional  navigation  and  then  switches  to  bang-bang 
when  the  missile  approaches  within  2km  from  the  target. 

A  summary  of  the  miss  distances  is  tabulated  below  for  the  two  scenarios  considering 
different  drag  models. 


Miss  Disi 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

0.4969 

18.6431 

0.2782 

3.0831 

VCPN 

0.5420 

18.1240 

0.0605 

2.8216 

Bang-bang 

0.3318 

0.5214 

Bang-PN 

MPOTTHHMl 

3.8627 

PN-bang 

0.5422 

0.5369 

0.4403 

0.3169 

Table  5a.  Comparison  of  Miss  Distances  (Tt  =  1.0s  and  Tm  =  0.5s)  without 
considering  drag  due  to  missile  turn-rate 


Miss  Disi 

tance  (m) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

0.1290 

19.8792  • 

0.1381 

4.5232 

VCPN 

0.0480 

19.0654 

0.1240 

3.9718 

Bang-bang 

0.5531 

0.0854 

Bang-PN 

7.2254 

PN-bang 

0.0512 

0.7763 

0.3128 

0.2759  ! 

Table  5b.  Comparison  of  Miss  Distances  (Tt  =  1.0s  and  Tm  =  0.5s)  with  additional 
drag  due  to  missile  turn-rate 


Bang-PN  Strategy 

We  observed  that  the  Bang-PN  strategy  produced  the  largest  miss  distances 
against  maneuvering  targets.  Switching  to  PN  guidance  at  the  terminal  stage  when  the 
target  maneuvers  is  not  effective.  In  fact,  it  is  when  the  missile  is  nearer  to  the  target, 
that  larger  accelerations  are  often  required  from  the  missile  to  achieve  target 
interception.  This  strategy  is  thus  not  effective. 
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PN-Bang  Strategy 


On  the  other  hand,  we  find  that  the  miss  distances  obtained  using  the  PN-Bang 
strategy  against  the  9-g  target  maneuvers  were  much  lower  than  PNG,  VCPN  and 
Bang-PN,  and  comparable  to  Bang-bang’s  results. 

A  comparison  of  the  amount  of  missile  control  effort  in  given  in  Table  6.  We 
observe  that  the  control  effort  required  by  the  PN-Bang  strategy  is  almost  consistently 
the  lowest. 


Missile  Control  Effort  (1 

ntegration  of  turn-rate) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

76.2543 

200.9011 

622.9764 

711.1141 

VCPN 

75.502 

199.255 

608.6091 

681.0754 

Bang-bang 

149.3850 

243.9412 

671.2259 

701.4282 

Bang_PN 

145.5924 

253.6009 

678.3973 

769.8148 

PN_Bang 

75.5047 

193.2745 

640.6799 

667.8471 

Table  6a.  Missile  Control  Effort  (Tt  =  1.0s  and  Tm  =  0.5s)  without  considering  drag 
due  to  missile  turn-rate 


Missile  Control  Effort  (1 

[ntegration  of  turn-rate) 

Scenario#! 

Scenario#2 

No  maneuver 

9-g  maneuver 

No  maneuver 

9-g  maneuver 

PNG 

76.2628 

196.771 

651.076 

736.5428 

VCPN 

75.4792 

196.0564 

609.9731 

698.861 

Bang-bang 

144.6337 

240.7891 

655.8412 

707.0795 

Bang_PN 

145.7526 

246.2364 

639.9804 

716.6012 

PN-Bang 

75.6336 

192.4796 

656.9838 

692.0159 

Table  6b.  Missile  Control  Effort  (Tt  =  1.0s  and  Tm  =  0.5s)  with  additional  drag  due  to 
missile  tum-rate 


Noise  in  LOS  rate 

Next,  we  examine  the  performance  of  the  PN-Bang  strategy  in  a  noisy  LOS  rate 
environment.  The  results  are  tabulated  in  Table  7. 
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Standard 

Deviation 

Miss  Distance  (m) 

Scenario#!  (9-g  maneuver) 

Scenario#2  (9-g  maneuver) 

Drag  w/o  turn 
rate  factor 

Drag  factoring 
turn-rate 

Drag  w/o  turn 
rate  factor 

Drag  factoring 
turn-rate 

(7=0  (noise  free) 

0.7763 

0.2759 

<7=0.001 

1.2190 

BWr^Tlff 

0.2998 

(7=0.005 

1.3813 

1.5121 

0.3274 

0.0432 

(7=0.01 

0.3762 

<7=0.02 

7.6481 

8.4838 

0.2061 

0.0473 

0=0.05 

40.2425 

41.3410 

21.5640 

6.3360 

Table  7.  PN-Bang:  Miss  distance  with  noisy  LOS  rate  (Tt  =  1.0s  and  Tm  =  0.5s) 


We  observe  that  the  performance  of  PN-Bang  strategy  is  tolerant  to  additive  white 


noise  in  the  LOS  rate  up  to  o=  0.01  rad/s. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  PROPORTIONAL  NAVIGATION 

PNG  is  known  to  be  optimal  against  a  non-maneuvering  target,  but  will  result  in  large 
miss-distances  against  evasive  target  maneuvers. 

2.  VELOCITY  COMPENSATED  PROPORTIONAL  GUIDANCE 

VCPN  was  found  to  offer  little  improvement  in  the  guidance  performance.  The 
conventional  PNG  has  a  closed  loop  update  on  the  closing  velocity,  Vc,  which  invariably 
provides  update  on  the  change  in  missile  velocity.  Hence  conventional  PNG  in  essence 
already  factors  a  correction  for  reduced  missile  velocity  into  the  guidance  acceleration 
command. 

For  the  case  of  a  non-maneuvering  target,  PNG  is  known  to  be  optimal,  hence  VCPN 
has  no  advantage.  However,  in  the  case  of  a  maneuvering  target,  the  additional  velocity 
compensation  term  happens  to  provide  a  small  additional  acceleration  at  the  final  moment  of 
target  interception  to  help  bring  about  a  slight  improvement  in  the  miss  distance  of  0.2-0.8m. 
This  is  still  not  good  enough  to  effectively  intercept  a  maneuvering  target. 

3.  BANG-BANG  GUIDANCE 

Bang-bang  guidance  produced  much  smaller  miss  distances  but  was  found  to  be 
sensitive  to  changes  in  LOS  rate,  especially  when  the  target  is  farther  away.  This  is  not 
desirable  in  the  initial  pursuit  of  the  target,  as  it  demands  more  control  effort  from  the  missile, 
which  swings  in  either  direction  of  the  LOS  due  to  the  full  20-g  swings  in  the  guidance 
commands.  This  also  makes  the  Bang-bang  strategy  vulnerable  to  noise  in  LOS  rate. 
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4.  COMBINED  PROPORTIONAL  NAVIGATION  WITH  BANG-BANG 

A  combined  strategy  incorporating  proportional  navigation  and  bang-bang  guidance 
was  investigated.  It  was  found  that  the  results  were  not  good  if  we  start  with  Bang-bang  and 
then  switch  to  PN  during  the  terminal  stage. 

The  better  strategy,  which  was  shown  to  be  effective  against  a  maneuvering  target,  was  PN- 
Bang,  which  started  the  target  pursuit  with  PN  guidance  and  then  switched  to  Bang-bang 
when  the  missile  was  nearer  to  the  target. 

This  strategy  worked  well,  since  the  change  in  LOS  rate  was  small  when  the  target 
was  farther  away  and  PN  guidance  was  adequate,  and  had  the  advantage  of  less  susceptibility 
to  noise.  When  the  missile  was  closer  to  the  target,  the  LOS  rate  change  was  expected  to  be 
larger,  especially  when  the  target  maneuvered,  and  Bang-bang  performed  well.  It  was  also 
observed  that  the  missile  control  effort  required  by  the  PN-Bang  strategy  was  similar  to 
conventional  proportional  navigation. 

B.  RECOMMENDATIONS 

This  strategy  is  thus  proposed  as  an  improvement  to  conventional  proportional 
navigation  guidance  against  an  evasive  target. 

This  study  explored  the  results  for  two  typical  engagement  scenarios.  Follow-on 
studies  can  be  carried  out  to  test  a  wider  engagement  envelope  for  this  guidance  strategy. 
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APPENDIX  A.  SIMULINK™  MODELS 


1.  Proportional  Navigation:  (Set  K=0  and  Guidance  Algorithm=PNG.m) 

2.  Velocity  Compensation  Proportional  Navigation:  (Set  K=1  and  Guidance 
Algorithm=PNG.m ) 

3.  Bang  Bang:  (Set  K=0  and  Guidance  Algorithm  =  Bangbang.m) 
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4.  Hybrid  Bang  Bang  with  VCPN: 

a)  Starts  with  Bang-bang,  switches  to  VCPN  at  R<2km  (Bangpn.m) 

b)  Starts  with  VCPN,  switches  to  Bang-bang  at  R<2km  (PNBang.m) 


Hybrid  Bang-bang  with  VCPN  (Vc=Rdot)  Implementation 


TgtTumRate  I 


MATLAB 

Function 


Target  Dynamics 


Integrator 


Random 


Number 


ThetaLdotnoisy 


GuidancePN 


|omegaout[^ — 
Missile  Turn  Rate 


Transfer  Fen 


MATLAB  . 

I —  <4—  Mux 

Function 

Velocity  Compensation 
Vmdot.an(ThetaML)  Mux5 


MATLAB 

Function 

Missile  Dynamics 


Integrator! 


Missile  Acc  X 


Absl  Integrator 


Integrated  Ms!  Accl 


Missile  Acc  Y 


Missile  Blocks: 


Missile  Dynamics: 


MatLab  Function  ‘MissileDynamics.m’ 
Missile  Integrator:  Initial  value  =  minit 


MissileOut:  Variable  Name=MissileOut 

Max  No.  of  Rows  =  OutputVec 


Guidance  Blocks: 


Matlab  Functions: 

a)  Pnvariables.m:  -  Computes  necessary  variables 

b)  GuidancePN:  -  Guidance  Algorithm 

c)  Velcomp.m:  ~  turnrate  compensation  for  missile 
velocity  due  to  drag 


Target  Blocks: 

T gtT urnrate  Step :  Step  time: 

i)  No  turn  =  tmax+1 

ii)  9g  (9x9. 81/Vt)  turn  =  intercept  time  -  2 
Target  Dynamics: 

MatLab  Function ‘TwoDdynamics.m’ 

Target  Integrator:  Initial  value  =  tinit 


TgtOut:  Variable  Name=MissileOut 

Max  No.  of  Rows  =  OutputVec 
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APPENDIX  B.  SOURCE  CODES 


The  list  of  Matlab  source  codes  are: 


S/N 

Filename 

Task 

1. 

Init.m 

Initializes  the  following: 

a)  missile  and  target  state  vectors 

b)  Missile  Drag  variables 

c)  Missile’s  initial  launch  velocity 

d)  Sampling  time  and  maximum  missile  flight  time 

2. 

Stop.m 

Stops  the  simulation  run  when  missile  passes  target. 
This  occurs  when  closing  velocity  is  +ve  (set  at 
Vc>=100m/s). 

3. 

Missiledynamics.m 

Computes  missile  state  vector  given  the  commanded 
turn  rate. 

It  considers  the  deceleration  in  the  missile’s  coasting 
speed  due  to  aerodynamic  drag.  The  drag  also  takes 
into  account  the  effect  of  turn  rate. 

4. 

T  argetdynamics .  m 

Computes  the  target  state  vector  given  the  commanded 
turn  rate.  No  drag  component  is  considered  as  it  is 
assumed  that  the  aircraft  has  sustained  thrust  to 
overcome  its  drag  components. 

5. 

PNvariables.m 

Computes  the  variables  for  Proportional  Navigation 
missile  guidance. 

Inputs:  Missile  and  Target  State  Vectors. 

(X-pos,  X-vel,  Y-pos,  Y-vel) 

Returns:  [ThetaLdot,R,ThetaML,Vm] 

6. 

Guidance  Algorithm: 

a)  PNG.m 

b)  Bangbang.m 

c)  BangPN.m 

d)  PNbang.m 

Computes  the  guidance  command  (msl  turn  rate) 
subject  to  20-g  limit. 

a)  Conventional  True  Prop.  Nav. 

b)  Bang-bang  guidance. 

c)  Starts  with  20-g  Bang-bang  and  switches  to  Prop 

Nav  with  N’=3  when  target  is  within  2km. 

d)  Starts  with  Prop  Nav  with  N’=3;  switches  to  20-g 
Bang-Bang  when  target  is  within  2km. 

7. 

Velcomp.m 

Computes  the  additional  missile  turn  rate  to 
compensate  for  reduced  speed  of  missile  due  to  Drag. 

8. 

Plotfigures.m 

Plots: 

Fig.  1:  Flight  trajectories  of  missile  and  target 

Fig.  2:  Miss  distance 

Fig.  3:  Missile  Speed  profile 

Fig.  4:  Missile  acceleration  profile 

Fig.  5:  Missile  acceleration  control  effort  (Integration 
of  acceleration  commands) 
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1.  Initialization 


%%  Init.m 

%  This  script  file  initializes  the  point  mass  missile  intercept  simulation 
global  Vcinit  beta 

tmax=ll;  %  Seconds  -  max.  total  simulation  time 

stime  =  0.001;  %  Sampling  time 

OutputVec= [20000,1, stime] ;  %  For  output  vector 


%  Computing  for  Missile  Drag 

h=2 000 ;  %  height  of  missile  during  cruise 

radius=0.15;  %  radius  of  missile  body 

Area=pi*radiusA2 ;  %  Cross-sectional  Area  of  missile 
CD=0.2;  %  Coefficient  of  Drag 

Mas s  =  3  00 ; 

if  h<9144  %  below  9144m 

rho=l . 22557 *exp ( -h/9144 ) ; 
else 

rho=l. 75228763 *dxp ( -h/6705 . 6) ; 

end 

beta=rho*CD*Area/ (2*Mass) ;  %Drag=beta*VmA2; 

%Vmdot= Drag /Mass  (deceleration  of  Vm) 

%  Initialize  Target  State  Vector 
%  tl=  scenario#!  and  t2=scenario#2 


Vt=500 ; 


tl=[  10000;  % 

-Vt ;  % 

1000;  % 

0];  % 

t2=  [  8000;  % 

-Vt*sin (pi/4) ;  % 

800;  % 

Vt*sin (pi/4) ] ;  % 


scenario=  input ( ' Enter 


x-position 

x-velocity 

y-position 

y-velocity 

x-position 

x-velocity 

y-position 

y-velocity 

inario  no.  "1"  or  "2"  :  '); 


if  scenario  ==  1 
tinit=tl; 
else 

if  scenario==2 
tinit=t2 ; 
else 

scenrario=  input ( 1  Please  re-enter  scenario  no.  "1"  or  "2"  :  ' ); 

end 

end 
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%  Initialize  Missile  State  Vector 
%  Missile  is  initially  twice  the  speed  of  target 


Vm=2  *Vt ; 


minit=[  0;  % 

Vm*cos (atan2 (tinit (3) , tinit (1) ) ) ;  % 

0;  % 

Vm*sin (atan2 ( tinit (3 ), tinit ( 1) ))] ;  % 


%  Calculate  initial  closing  velocity  for  Approx#2 


x-position 

x-velocity 

y-position 

y-velocity 


LOSinit=atan2 (tinit (3) , tinit ( 1 ) ) ; 
Vminit=sqrt (minit (4) /v2+minit (2)A2); 
Vtinit=sqrt ( tinit ( 4 ) ^2+ tinit (2)^2); 
ThetaMinit=atan2 (minit (4) , minit (2) ) ; 
ThetaTinit=atan2 (tinit (4 ) , tinit (2 )  )  ; 


%  Initial  LOS  Target  to  Missile 
%  Initial  Missile  velocity 
%  Initial  Target  velocity 
%  Initial  Missile  heading 
%  Initial  Target  heading 


Vcinit=Vtinit*cos (ThetaTinit-LOSinit ) -Vminit*cos (ThetaMinit-LOSinit ) 


2.  Stopping  Simulation 

%  Stop.m 

%  Stop  Simulation  if  Vc--fve  ...  ie  missile  pass  target. 

function  [stop]  =  stop(u) 

%  Checks  if  Vc  is  closing  or  not 
if  u ( 1 ) >=100 
Stop=l ; 
else 

stop=0 ; 

end 
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3.  Missile  Dynamics 


%  Missiledynamics.ru  :  Mat  lab  Function  used  in  Simulink  model  "  PNBang .  mdl 
%  Computes  the  2 -Dimension  Missile  dynamics . 

%  Inputs:  State  Vector  (X-position,  X-velocity,  Y-position,  Y-velocity) 

%  Returns:  Missile  Turn  Rate 

function  [xdot]  =  Missiledynamics (u) 

global  beta 

w=u(l) ;  %  turnrate 

x=u(2:5);  %  state  vector  :  x(t) 

Vm=sqrt ( (x(2) ^2+x(4) A2) ) ;  %  Missile  velcoity,  Vm; 


A= [ 0 ,  l,  0,  0; 
0,  0,  0 ,  -w; 
0,  0,  0,  1; 
0,  wf  0,  0]; 


%  Compute  Drag 

%  Nominal  frontal  Drag  due  to  aerodynamic  surface 
Drag0= [ 0 ;  x (2 ) *beta*Vm;  0;  x ( 4 ) *beta*Vm] ; 

%  Assuming  Drag_max  is  about  SxDragO  at  maximum  turnrate  (i.e.  w_max ) 

w_max=20*9 . 81/Vm;  %  Max.  turn  rate  =  20g/s 

k=4/  (w_max^2 )  ; 

Drag=Drag0* (l+k*w^2 ) ;  %  Drag  is  a  function  of  turnrate. 

%  Output  state  vector  :  xdot=A(v^)  .  x  -  Drag 
xdot=A*x-Drag; 


4.  Target  Dynamics 

%  Target dynamics .m  :  Matlab  Function  used  in  Simulink  model  " PNBang .mdl ” 
%  Computes  the  2-Dimension  Target  dynamics. 

%  Inputs:  State  Vector  (X-position,  X-velocity,  Y-position,  Y-velocity) 

%  Returns:  Target  Turn  Rate 

function  [xdot]  =  PNdynamics (u) 

w=u(l) ;  %  turnrate 

x=u (2:5) ;  %  state  vector  :  x(t) 

%  Output  state  vector  :  xdot=A(v;)  .x 

A=  [  0 ,  1,  0,  0; 

0,  0,  0,-w; 

0,  0,  0,  1; 

0,  w,  0,  0]  ; 

xdot=A*X; 
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5.  Proportional  Navigation  Variables 


%  PNva r i ab 1 e s . m 

%  Mat  lab  Function  used  in  Simulink  model  11  PNBang .  mdl  " 

%  Computes  the  variables  for  Proportional  Navigation  missile  guidance. 

%  Inputs:  Missile  and  Target  State  Vectors . (X-pos,  X-vel,  Y-pos,  Y-vel) 
%  Returns :  [ThetaLdot , R, ThetaML, Vm] 

function  [output]  =  PNvariables (u) 

%  Compute  the  LOS  (missile- target )  turnrate 
TgtOut=u (1:4) ; 

MissileOut=u (5 : 8) ; 

xt=TgtOut ( 1) ;  %  Target  x-position 

Vxt=TgtOut (2 ) ;  %  Target  x-direction  velocity 
yt=TgtOut (3 ) ;  %  Target  y-position 

Vyt=TgtOut (4) ;  %  Target  y-direction  velocity 

xm=MissileOut (1) ;  %  Missile  x-position 
Vxm=MissileOut (2) ; %  Missile  x-direction  velocity 
ym=MissileOut (3 ) ;  %  Missile  y-position 
Vym=MissileOut (4) ; %  Missile  y-direction  velocity 

%  Compute  ThetaL:  LOS  angle  between  missile  and  target 
ThetaL=atan2 ( (yt-ym) , (xt-xm) ) ; 

%  Compute  The taM:  Missile  heading 
ThetaM=atan2 (Vym,Vxm) ; 

%  Difference  of  angle  between  Missile  Heading  to  LOS  (Msl-Tgt) 

ThetaML  =  ThetaM-ThetaL; 

%  Compute  R:  Range  of  target  to  missile 
Rsquare= (xt-xm) A2+ (yt-ym) A2 ; 

R=sqrt (Rsquare) ; 

%  Missile  velocity 
Vm=sqrt  (VxmA2+VymA2 )  ; 

%  Compute  The t ado t 

ThetaLdot=  (  (xt-xm)  *  (Vyt-Vym)  -  (yt-ym)  *  (Vxt-Vxm)  )  /Rsquare; 

%  PN  variables  output 

output = [ThetaLdot ; R; ThetaML; Vm] ; 
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6.  Guidance  Algorithm 


a)  PNG.m:  Conventional  True  Proportional  Navigation 

%  PNG.m 

%  Matlab  Function  used  in  Simulink  file  "PNguidance" 

%  Computes  the  guidance  command  {missile  turn  rate)  for  True  Prop  Nav 
%  Inputs:' 

%  u(l)  =  ThetaLdot  (rate  of  change  of  LOS) 

%  u  (2 )  =  Raot=Vc  (closing  velocity) 

%  u(3)=  ThetaML  (missile  leading  angle) 

%  u(4)  =  Vm  (missile  velocity) 

%  u(5)=  R  (Range  of  target  from  msl) 

%  Returns:  Missile  Turn  Rate  based  on  PN 

function  [omegaout]  =  PNG(u) 

%  Compute  the  missile  turnrate 

omegaoutl=-u (2 ) *u  (1) ;  %  Vc*ThetaLdot/Vm* cos (ThetaML) 

%  PropNav:  N' *Vc* ThetaLdot /Vm* cos (ThetaML) 
omega=3*omegaoutl/ (u (4) *cos (u (3)  )  )  ; 

%  20-g  limit 

g„limit=20*9 . 81*sign (omegaoutl) /u(4) ; 
if  abs (omega)  >=  abs(g_limit) 
omegaout=g_limit ; 
else 

omegaout =omega ; 

end 


b)  Bangbang.m:  20-g  Bang  bang  guidance 

_%  Bangbang - m 

%  Matlab  Function  used  in  Simulink  file  "BangBang" 

%  Computes  the  guidance  command  (missile  turn  rate)  for  PN  Bang-bang 
missile  guidance . 

%  Inputs : 

%  u ( 1 ) =  ThetaLdot  (rate  of  change  of  LOS) 

%  u ( 2 ) =  Rdot=Vc  (closing  velocity) 

%  u{3)=  ThetaML  (missile  leading  angle) 

%  u(4)=  Vm  (missile  velocity) 

%  Returns:  Missile  Turn  Rate  based  on  PN 

function  [omegaout]  =  Bangbang (u) 

%  Compute  the  missile  turnrate 

omegaout l=-u (2) *u(l) ;  %  Vc -ThetaLdot /Vm* cos (ThetaML) 

omegaout=20*9 . 81*sign (omegaoutl ) /u (4) ; 
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c)  BangPN.m:  Starts  with  Bang-bang;  switches  to  Prop  Nav  at  Range<2km 


%  Bang pn . m 

%  Matlab  Function  used  in  Simulink  file  "Hybrid  BangPN" 

%  Inputs: 

%  u(l)  =  ThetaLdot  (rate  of  change  of  LOS) 

%  u(2)=  Rdot=Vc  {closing  velocity) 

%  u ( 3 ) =  ThetaML  (missile  leading  angle) 

%  u(4)=  Vm  (missile  velocity) 

%  u(5)=  R  (Range  of  target  from  missile) 

%  u(6)=  Vmdot sin (ThetaML) /cos (ThetaML)  =  Velocity  compensation 
%  Returns:  Missile  Turn  Rate  based  on  PN 

function  [omegaout]  =  Bangpn(u) 

%  Compute  the  missile  turnrate 
%  For  TPN  Implementation 

omegaout l=-u  (2 )  *u  ( 1)  ;  %  Vc*ThetaLdot/Vm* cos  (ThetaML.) 

if  u ( 5 ) <2000 

omega=3 * omegaout 1/ (u (4 ) *cos (u ( 3 ) ) ) -u ( 6 ) ;  %  PropNav: 

N' *Vc*ThetaLdot/Vm*cos (ThetaML) 

%  20-g  limit 

g_limit=20*9 . 81*sign (omegaoutl) /u (4) ; 
if  abs (omega)  >=  abs(g_limit) 
omegaout = g_l imi t ; 
else 

omegaout =omega ; 

end 

else 

omegaout=20*9 . 81*sign (omegaoutl ) /u (4 ) ;  %  bang  bang 


d)  PNBang.m:  Starts  with  Prop  Nav  and  switches  to  Bang-bang  at  Range  <2km 

%  PNBang.m 

%  Matlab  Function  used  in  Simulink  file  "Hybrid  PNBang" 

%  Inputs: 

%  u ( 1 ) =  ThetaLdot  (rate  of  change  of  LOS) 

%  u ( 2 ) =  Rdot=Vc  (closing  velocity) 

%  u(3)=  ThetaML  (missile  leading  angle) 

%  u(4)=  Vm  (missile  velocity) 

%  u(5)=  R  (Range  of  target  from  missile) 

%  u(6)=  Vmdot sin (ThetaML) /cos (ThetaML)  =  Velocity  compensation 
%  Returns:  Missile  Turn  Rate  based  on  PN 

function  [omegaout]  =  PNbang(u) 

%  Compute  the  missile  turnrate 
%  For  TPN  Implementation 

omegaout l  =  -u (2) *u(l)  ;  %  Vc*ThetaLdot/Vm*cos (ThetaML) 

if  u ( 5 ) >2 000 

omega=3 * omegaoutl/ (u (4) *cos (u (3 ) ) ) -u ( 6) ;  %  PropNav: 

N' *Vc*ThetaLdot/Vm*cos (ThetaML) 
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%  20-g  limit 

g_limit=20*9 . 81*sign ( omegaoutl ) /u(4)  ; 
if  abs (omega)  >=  abs(g_limit) 
omegaout=g_limit ; 
else 

omegaout = omega ; 

end 

else 

omegaout=20*9 . 81*sign (omegaoutl) /u (4 ) ;  %  bang  bang 


7.  Velocity  Compensation 


%  Velcomp.m 

%  Matlab  Function  used  in  PN  Bang-Bang  guidance 

%  Computes  the  additional  missile  turn  rate  to  compensate  for  reduced  Vrn 
%  Inputs:  u(l)=  Vm  (missile  velocity) 

%  u(2)=  Vmdot  (missile  deceleration) 

%  u(3)~  ThetaML  (missile  leading  angle) 

%  Returns:  Missile  Turn  Rate  Compensation 

%  [Vmdot* sin (ThetaML) /Vm*cos (ThetaML) ] 


function  [comp]  =  Velcomp(u) 

%  Compute  the  missile  turnrate  compensation  due  to  reduction  in  Vm 
comp=u(2) *tan(u(3) ) /u(l) ;  %  Vmdot* sin (ThetaML) /Vm*cos (ThetaML) 


8.  Plotting 

%%  Plotf igures  .m 


% 

To  plot 

: 

'O 

Fig. 

I : 

Flight  trajectories  of  missile  and  target 

% 

Fig. 

2: 

Miss  distance 

% 

Fig. 

3: 

Missile  Speed  profile 

% 

Fig. 

4  : 

Missile  acceleration  profile 

% 

Fig. 

5: 

Missile  acceleration  control  effort  (Integration  of  acceleration 

commands ) 

Missilepos= [MissileOut (  : ,  1)  , MissileOut ( :  ,  3  )  ]  ; 
Tgtpos= [TgtOut (  : , 1)  , TgtOut (  :  ,  3 )  ]  ; 
time=MissileOut {  : ,  5 )  ; 

Mslaccx=accx/9 .  81  ; 

Mslaccy=accy/9 .81; 

%  Compute : 

%  a)  Miss  Distance 
%  b)  Vm 

%  c)  Missile  Accelerations 

missdis= [ ] ; 

Vm=  [  ]  ; 
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Mslacc=  [  ]  ; 


for  i=l :max ( size (MissileOut ) ) , 

%  a)  Miss  Distance 

mi ss=norm (Missilepos (i, : ) -Tgtpos (i, : ) ) ; 
missdis= [missdis ,miss]  ; 

%  b)  Vm 

Vml= (MissileOut (i, 2) ^2+MissileOut ( i # 4 ) ^2 ) ^0 . 5; 
Vm= [Vm, Vml ] ; 


%  c)  Missile  Accelerations 
Am=sqrt (Mslaccx(i) A2+Mslaccy (i) ^2) ; 
ThetaM=atan2 (MissileOut (i ,  4 ) , MissileOut (1,2)); 
Theta=atan2 (Mslaccy (i) ,Mslaccx(i) ) ; 
acc=Am*sin (Theta-ThetaM) ; 

Mslacc= [Mslacc , acc]  ; 


end 

mindis=rain  (missdis)  %  Minimum  Miss  Distance  (i.e.  Closest  Pt  of  Approach) 
index=f ind (missdis==mindis ) ;  %  Index  where  min.  distanc  occurs 

intercept=sqrt (Tgtpos ( index, 1 ) ^2+Tgtpos ( index, 2)^2);  %  Intercept  point 
Controlef fort=IntgAcc ( 1 : index) . *Vm( 1 : index) %  Acc  Control  Effort 


figure (1) 
elf 

plot (Tgtpos ( : , 1) , Tgtpos ( : , 2) , 'r' ) 
hold  on 

plot (Missilepos ( : ,1) ,Missilepos( : ,2) ) 

plot (Tgtpos (index, 1) , Tgtpos (index, 2) , 'x' ) 

plot (Missilepos (index, 1) , Missilepos ( index, 2 ) , 'o' ) 

title ('Rear  Cross-Range  Scenario  Geometry') 

xlabel ( 'x-coordinates  (meters) ' ) 

ylabel ( 'y-coordinates  (meters) ' ) 

gtext ( 'Missile' ) 

gtext ( 'Target ' ) 

gtext ([' Intercept  range  =', num2str ( intercept )] ) 

figure (2) 
cl  f 

plot ( time, missdis) 

title ('Rear  cross-Range  Scenario:  Miss  Distance') 

xlabel ('time  (sec)') 

ylabel ( 'Miss  Disttance  (meters) ' ) 

gtext ([ 'Miss-distance (m) =' , num2str (mindis ) , '  Time  of  Intercept 

(s) = ' , num2str ( (index-1) *stime) ] ) 

figure (3 ) 
cl  f 

plot  (time,  Vm) 

title ('Rear  cross-Range  Scenario:  Missile  Velocity') 
xlabel ( ' time  (seconds)  '  ) 

ylabel ( 'Missile  vector  velocity,  Vm  (meters/sec)') 
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grid 


figure (4) 
elf 

plot (time (1 : index) , Mslacc (1 : index) ) 
title ( 'Missile  Absolute  Acceleration') 
xlabel ( ' time  ( sec ) ' ) 
ylabel ( 'Acceleration  in  g') 

%gtext {[ 'Missile  Initial  Accelertion  =  ', num2str (Mslacc ( 1 ))] ) 

gtext ( [ 'Missile  Acceleration  Effort=' , num2str (Controlef fort ( index)  )  ]  ) 

grid 


figure (5) 
elf 

plot  (timed:  index)  #  Mslaccx  (1 :  index)  ,  ,  time  (1 :  index)  ,  Mslaccy  ( 1 :  index)  , 

legend ( 'accx' , 'accy' ) 

grid 

title ( 'Missile  Acceleration  in  x-y  directions') 
xlabel ( ' time  ( sec ) ' ) 
ylabel ( 'Acceleration  in  g') 

figure (6) 
elf 

plot (time ( 1 : index) , Controlef fort , ' - ' ) 
title ( 'Integrated  Missile  Control  Effort') 
xlabel ( ' time  ( sec ) ' ) 
ylabel ( 'Acceleration  in  g') 

gtext ( [ 'Missile  Acceleration  Effort^' , num2str (Controlef fort (index) )  ]  ) 
hold  on 
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APPENDIX  C.  SIMULATION  RESULTS 


This  appendix  contains  the  following  plots  for  the  different  guidance  strategies: 

Fig.  1 :  Flight  trajectories  of  missile  and  target 

Fig.  2:  Miss  distance 

Fig.  3:  Missile  speed  profile 

Fig.  4:  Missile  acceleration  profile  and  control  effort 

The  guidance  strategies  are: 

1 .  PN  with  different  turn-rate  time  constants 

2.  VCPN 

3.  Bang-bang 

4.  Bang-bang  with  PN:  Starts  with  Bang-bang  switching  to  PN  when  R<2km 

5.  PN  with  Bang-bang:  Starts  with  PN  switching  to  Bang-bang  when  R<2km 
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Fig.  3  Missile  Velocity 
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Scenario  #2 
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Scenario  #2  (with  additional  drag  due  to  turn-rate) 
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Fig.  4  Missile  Acceleration  Profile 
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Fig.  4  Missile  Acceleration  Profile 
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ThetaLdot:  Effects  of  noise 
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Missile  vector  velocity,  Vm  (meters/sec) 


Missile  Velocity 


APPENDIX  D.  MATLAB®  INFORMATION 


MATLAB®  and  SIMULINK™  is  a  product  of  Math  Works,  Inc.,  24  Prime  Way, 


Natick,  Mass.  01760. 

MATLAB®  version  5.3  and  SIMULINK™  were  used  throughout  this  study. 
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